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ABSTRACT 


In recent years, the concept of using the rudder for roll stabilization of a ship in 
a seaway has been investigated with good results. Such designs have been used to 
solve the roll problem of a ship on steady course while compensating for yaw 
perturbations. 

To complete a review of the possible design strategies to meet the specifications 
for this model, the effects of the feedback gains on rolling and yawing are studied in 
detail. Roll angle and roll rate feedback are used to control the rudder. 

Roll stabilization with the rudder in various sinusoidal sea states is studied by 
simulation on the IBM digital computer. The model used is based on the data 
obtained from a typical naval ship. The Root Locus method is used to design the 
feedback gains. The computer simulation programs are written in Digital Simulation 
Language (DSL/VS), are plotted as data in DISSPLA and include the effects of rudder 


servo nonlinearities, which seriously restrict the ability of the rudder to reduce roll. 
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I. INTRODUCTION 


There is a need to reduce roll on various types of ships in order to increase the 
comfort of the passengers, to provide safe operation of container ships in bad weather 
and to provide stable platforms for helicopter operations. 

Passive or active tanks, or stabilization fins are most commonly used at present 
to control roll. But both have disadvantages: tanks require a lot of space, fins 
introduce a considerable drag, are expensive, and also require space for the hydraulic 
actuators. 

Recently there has been a considerable amount of interest in Roll Stabilization 
systems which use the rudder (RRS), since rudder motions not only affect a ship's 
heading but influence the roll motions as well. This resulted in the successful 
installation of a rudder roll stabilization system on some ships of the United States 
Coast Guard Cutter HAMILTON class (378 foot) as reported by references (1) and 
cy 

This thesis describes a simple mathematical model for the transfer between the 
rudder angle and the two outputs : Rate of turn (heading rate) and roll angle. 
Computer simulations demonstrate results of a design procedure for a combined 
controller for roll stabilization with the rudder. Models are also provided for the 
disturbances and the steering machine. 

This thesis is organized as follows : Chapter II gives mathematical models of the 
ship, the disturbances and the steering machine. Chapter III deals with the controller 
design. Chapter IV gives results of computer simulation and Chapter V summarizes the 


results giving conclusions and recommendations for further research. 
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Il. MATHEMATICAL MODELING 


The motions of a ship in waves depend on 


e The dynamics of the ship 
e The disturbances and 
e The controller output, which is influenced by the steering machine. 


A. THE SHIP’S DYNAMICS 
In van Amerongen and Van Cappelle, 1981 [Ref. 3: p. F2 1-3], the basic 
mathematical model was derived which describes the dynamics between the rudder as 


input signal and yaw and roll as output signals. The results will be summarized below. 


The basic equation are : 


Y = m(v-ur). (aly 
K=1, ©. (2.2) 
Ne= | i (253) 
where 
m is the mass of the ship, including the mass of the 


displaced water. 


[he U. + a Ww, : Linear velocity vector with components along x, y, 


and Z-axis. 
u 1s rate of surging. 
Vv is rate of swaying. 
Ww is rate of heaving. 
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82=p, tqitty : Vector angular velocity with components about 


X, y, and z-axis respectively. 


p is rate of roll (= @). 
q IS tatenomsOMCn)( = Q). 
is rate of yaw (= ‘P). 


k= X ty FZ: Vector force acting on the ship. 


X is hydrodynamic force along x-axis. 


= 


is hydrodynamic force along y-axis. 


Z is hydrodynamic force along Z-axis. 


M=K,+ M, +N, : Vector moment acting on the ship. 


K is rolling moment about x-axis. 
M is pitching moment about y-axis. 
N is yawing moment about z-axis. 


Angular moment=I, p+ i q;+ B ry. 


te is mass moment of inertia about x-axis. 
i is mass moment of inertia about y-axis. 
I is mass moment of inertia about z-axIs. 


The coordinate system where the above variables are defined 1s shown in Figure 2.1. 
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Figure 2.1 Coordinate System of the Ship. 


IS 


The left-hand sides of equation (2.1) - (2.3) can be expanded into a Taylor series 


in which the higher - order terms are neglected : 


econ Yp® + Y50 (2.4) 
K =Kyv+Kr+ Kg® + Kg® + K56 (2.5) 
N=Nv+Nr+ N@o® + N56 (2.6) 


Laplace transform of equations (2.1) - (2.6), substitution of equations (2.1) and 
(2.4) in (2.5) and (2.6) and substitution of equations (2.2) and (2.3) in (2.5) and (2.6) 
yield, under the assumption that v is small : 


» 


~ 


@M(s)= (Ks 0(s) - K. 1(s)}. (2.7) 


s* + 2fm sto? 
I 
S) En iNs 0(s) = No @M(s)}. (2.8) 


These equations can be illustrated with the block diagrams of Figure 2.2 and 
RACuUTe 2a, 

The model used in the subsequent analysis was “derived” as follows. First it was 
assumed that the influence of the rudder on heading rate would not be strongly 
influenced by roll angle. The conclusion of this assumption is that a heading rate 
model could be established using only rudder angle as input. The second assumption 
was that linear models would be sufficient to identify the major features of the 
dynamics that would be required. With the further assumptions that the roll dynamics 
could be modeled as a second order system with inputs of rudder and heading rate 
only, the trial model looks like equation (2.7) and (2.8). 
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Figure 2.2. Block Diagram for a and Control System 
(With Compensator). 


where 


o 


is rolling angle. 


= 


is heading rate. 


oO 


is rudder angle. 


Zz 
oO 


Is a constant relating rudder angle to heading rate. 


an 
oO 
a 


a constant relating rudder angle to rolling. 


Fe 
n 


a constant relating heading rate to rolling. 


No is a constant relating rolling rate to heading rate. 


a 
past 


is the time constant of the heading rate system. 
Si is the Laplace transform variable (= d/dt). 
C is the damping ratio of the roll dynamics. 
QO, is the natural frequency of the roll dynamics. 
The parameter values of this model have been estimated from full-scale trials 


with a naval ship. Table t gives some results (Ref. 4: p. 44]. Note that Nq@ has been 
assumed to be zero. 
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Figure 2.3. Block Diagram of the Dynamics between Rudder and 
Yaw and Roll (Without Compensator). 
TAB GE 
PARAMETERS FOR A NAVAL SHEIIP 


| Fat 13 knots | at 17.5 Kneis at 20.5 knots 





(Note : All the above coefficients assume that the basic unit of time is the second. The 


basic unit of angular displacement can be either radians or degrees without altering the 


coefficients. However, for convenience, in the remainder of this thesis all angles will be 
in degrees.) 

It is interesting to compare the coefficients obtained for the naval ship in Table 1 
[Ref. 4: p. 44], with those of a ship for which a Rudder-Roll-Stabilization has 
apparently been successfully implemented. Ihe comparison data contained in Table 2 
was obtained for the CONFIDENCE from John R. Ware [Ref. 1: p.11]. There is a 
remarkable similarity in most of the parameters, there is a slight difference in the 
coefficient relating rudder angle to turn rate. Clearly, the CONFIDENCE will turn 


much more rapidly, with a consequent reversal of the rudder induced roll. 


TABEE=2 
PARAMETERS FOR CONFIDENCE MODEL 


ee 
80 








B. THE DISTURBANCES 

Waves are the most important disturbance with respect to roll. They can be 
described by means of a frequency spectrum, for instance the Bretschneider Spectrum 
[Ref. 5], or the Pierson-Moskowitz Spectrum [Ref. 6]. However, these models of the 
sea state are not conveniently used for the design of a roll stabilizing compensator. In 
this thesis the wave disturbance used is a single frequency sinusoid, with @ = 0.5 
rad/sec which is the natural frequency of the ship in roll. This disturbance is used to 
cause roll motion, and the compensator is adjusted to provide maximum reduction in 
roll amplitude. 

Two methods are used to verify the behavior of the system over the range of 


frequencies in a normal sea state. These were 
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(1) A number of single frequency sinusoids, of different frequencies and properly 
adjusted amplitudes were used individually to determine the spectral response 
as in Figures 4.21 through 4.25. 

(2) To approximate a real sea, eight sine waves of appropriate frequencies and 
amplitudes were added and the resulting wave used to check the roll 


suppression, as shown in Figures 4.3 through 4.12. 


C. THE STEERING MACHINE 

In order to reduce the roll motion of a ship using its rudder, the steering system 
must meet certain economic, hydrodynamic and machinery criteria. These criteria can 
be met by using as much of the existing expensive steering system machinery as 
possible, thus voiding most of the capital and maintenance costs associated with 
conventional roll stabilizers. 

When the rudder is going to be used for reduction of roll motions it should be 
able to follow frequencies near the natural roll frequency, @,, without a noticeable 
phase lag. The steering machine used can be described by the simplified block diagram 


of Figure 2.4a, and autopilot can be described by the block diagram of Figure 2.4b. 





Figure 2.4a Bo ee of the Steering Machine 
KRSER = 10.0. 


A positioning system is used to position the rudder for small signal operation. 
Its bandwidth should be much larger than that of the ship. In order to prevent a 


phase lag, the angle of rotation and the angle rate have limits. 
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The rudder angle linut 1s determined either by the mechanical construction of the 
steering machine or is set by the autopilot. The rudder rate limit is determined by the 
construction of the steering machine and by the number of hydraulic pumps which are 
in Operation. 

In the design of the controller care must be taken to ensure that the derivative of 
the output signal of the controller is less than the maximum rudder rate itn order to 


prevent phase lag. 





= 
FQ 
ae 
= 
O 
ie 
ms 
Gi 


Block Diagram of the Autopilot 
0000.0, M1 = 0.025, @2= 1000.0. 
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Hf, CONTROLLER DESIGN 


A. INTRODUCTION 

The primary objective of this chapter is to present procedures for the design and 
compensation of the wave input to roll output control system. Compensation is the 
adjustment of the system in order to limit the maximum roll magnitude. 

The approaches to control system analysis and compensation used in this thesis 
are the Root-locus approach and the frequency-response (Bode) approach. The system 
was tested in 2 phases. In the first phase, the uncompensated system 1s tested by 
adjusting the frequency of the wave input. The wave function was generated as a 
sinusoidal input. Since both roll angle and roll rate are easily measured on ship, it was 
proposed that these be fed back to provide roll damping. 

The system which may be designed by a trial-and-error approach 1s checked to 


see Whether the designed system satisfies a desired maximum roll magnitude. 


B. ROOT-LOCUS APPROACH TO CONTROL SYSTEM DESIGN 
The Root-Locus method is a graphical method for determining the locations of 
all closed loop poles from knowledge of the location of the open loop poles and zeros 


as some parameter (usually the gain) 1s varied from Zero to infinity. 


G 


S 


p 
UJ x , 7 2. 
KN S“+2¢w, Stw_? 


— 


Gr Ge 


KRSER 


Figure 3.1 Block Diagram Manipulation for Root Locus Method 
ave Input to Roll Output). 
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In Figure 3.1 


on = 0.59 
C = 0.18 
KRSER = 10.0 


In Figure 3.1, the transfer function of the roll compensator, G_., should be 
G.(s)=K,*S+K, =K,(s+ K,/K,). (ea!) 


In order to solve this problem, the Root-locus method was used to check the 
feasibility of this proposal. First, the ship’s dynamics are expressed as a rolling transfer 
function. That is, the feedback term is neglected. 


In Figure 3.1, the roll to wave transfer function, G,, is 


. o7 0.3481 ce 
S ee ee : 
vst + 2G st+a,*  — s*+0.2124s + 0.3481 


Also, we can get the Root Locus for the roll transfer function. The result is 
eurero 2. 

The rudder servo is very fast compared to the natural frequency of the ship’s 
rolling mode, and the transfer function of the rudder servo, G., may be approximated 


as 


KRSER 10 
23) (253) 
; S+KRSER S+10 
The loop transfer function, G, re 1S 
3.481(K,s+K,) 
G,.(s)= a (3.4) 


(s+ 10)(s*+0.2124s + 0.3481) 


Note that the open loop has three poles and one zero. The poles are at fixed 
locations, and the location of the zero can be chosen as desired. This leads to a Root- 
Locus as shown on Figure 3.3, where the location of the zero has been chosen 


arbitrarily to illustrate the behavior of the roots. The actual location of this zero is 
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determined by the coefficients K, and K,. Root locations are changed by setting the 
gain K,. Proper choice of K, and K, determines the ¢ and @, of the roots and the 


damping of these roots 1s, of course, the damping of the ship roll motion. 


The closed loop transfer function, G,,, is 


0.348 1(s+ 10) 


G (8) = = D a 
s> + 10.21245 + (2.4721 + 3.481K,)s+3.481(1+ K,) 


(3.5) 


From equation 3.5: 

The zero at S=-10 should have little or no effect, since it is outside the 
bandwidth. Therefore the response depends on the roots of the denominator and on 
the. Ceca 


From equation 3.4: 

When we feedback roll + roll rate, we change the location of the zero. So the 
range of the zero is from “0” to “10.0” in calculation. For a small zero the root locus 
will have a loop near the origin, will go to the real axis, then out to the asymptote at 
about -5.0. When the zero is out near -1.0, the Root-Locus from the complex poles 
will not go to the real axis, but will stay complex as it goes to the asymptote. 

The best we can do is use the rudder to improve the damping of the roll. Then 


small zeros are not useful because we do not want to change @.. 
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Figure 3.2 Root Locus for Wave to Roll 
(Without compensator). 
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Figure 3.3b Root Locus for Stabilized Model (Main Part) 
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CF PRE OREN C TGR Sueno (BODE) APPROACH TO CONTROL SYSTEM 


In dealing with the problem of compensating control systems via frequency- 
domain techniques, we control the transient-response behavior in terms of such 
frequency-domain specifications as phase margin, gain margin, resonant peak value, 
and bandwidth. Design in the frequency domain is indirect because the system is 
designed to satisfy these frequency-domain specifications rather than time-domain 


specifications. 


Gs 


1 
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Figure 3.4 Block Diagram Manipulation for Bode Method 
(Wave Input to Roll Output). 
In Figure 3.4 
OY ae We, 
C = 0.18 
KRSER = 10.0 
From the block diagram on Figure 3.4, the closed loop contains a rudder servo, 
and a compensator block to stabilize the system and provide damping of system 


response to the wave input. The loop transfer function, G,,, is 


K {((K4/K)S+ 1} 


G,(9=—— St 
on” (O.1S + 1)(2.872782 + 0.61025 + 1) 


(3.6) 
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If K, > 1, then Bode gain increases and the gain crossover changes accordingly. 


The transfer function, equation 3.6, for G, oS) can be written as 


S(K,/K,) +1 l 


G 4 <<... lL 
(=ER ag (S2/0.3481) + (0.2124/0.3481)S + | 


}] (3.7) 


That is, the feedback (G,) and rudder servo (G_) transfer functions can be 
combined to have the algebraic form of a lead filter. Also, we can get the frequency 
responses (Bode diagram). The results are Figure 3.5 to Figure 3.7. Figure 3.5 shows 
the uncompensated roll to wave frequency response and Figure 3.6 and Figure 3.7 
show the roll to wave frequency responses (open loop and closed loop) with 
compensator (K, = 4.0, kK, = 4.0). From this, it was observed that the ship has a 
natural frequency of in roll 0.59 radians per second which corresponds to a period of 
10.65 seconds. 

Then the effect (seen on the Bode diagram in Figure 3.6) of any gain, K,, is to 
raise or lower the |G| curve. Raising the |G| moves the gain crossover to higher 
frequency, regardless of the location of the zero of the filter. 

The purpose of the filter zero is to increase damping of the roll frequency. To do 
this we want to introduce positive phase at the gain crossover. If we put the zero at a 
frequency lower than the gain crossover (caused by the gain K,), we increase the 
bandwidth which 1s undesirable. 

However, if we place the zero at a higher frequency, we increase the damping 
(phase margin) by a reasonable amount without appreciable increase in bandwidth. 


These results are essentially the same as shown by Root-Locus analysis. 


In Figure 3.6, the pole of the rudder servo becomes a Zero of the closed loop, and 
the effect of compensation is to move the real root and complex roots to new 
locations, presumably with larger negative real parts for the complex roots. If the 
increase in gain has been small, on the root locus the roots have not moved very far. 

Therefore the real effect of the compensator has been to move the complex roots. 
The effect can be seen on the closed loop Bode diagram. 

The above implies that roll would be reduced at low frequency, though not very 


much, but would be reduced substantially in the range of frequencies near resonance. 
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Figure 3.5 Frequency Response for Wave to Roll 
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IV. THE RESULTS OF COMPUTER SIMULATION 


A. SIMULATION TEST STUDIES 

This chapter will give a few results of the computer simulations. During the 
computer simulation, the parameter values used in this thesis are taken from in 
Table 1 with ship speed at 17.5 knots. The autopilot gain constant and rudder servo 
gain constant are given in Figure 2.4. 

The wave disturbance was generated as a sinusoidal input which produced a roll 
angle of + 20 degrees in the open system at the resonant frequency which 1s 0.5 
rad/sec. This was used for analysis and design. A complex wave consisting of a sum 
of sinusoids of different frequencies and amplitudes (see program, Appendixces E and 
F) was used to test the chosen designs. The test waves are shown in Figure 4.1. These 
input waves were used to drive the ‘open system’, the system with ideal rudder, and the 
system with real rudder. 

Figure 4.2 used the single sine wave to test the system with ideal rudder. 
Response of the compensated system is compared with that of the uncompensated 
system. In the ideal system with compensator in this reduces the maximum roll angle 
to approximately + 9 degrees. So the roll reduction is 55 percent. The numerical 
results are given in table 5 in Appendix I. System response to the complex wave 1s 
given in Figure 4.3. 

It appears that a considerable reduction with ideal rudder (approximately 71 
percent: Z = 2.0, K, = 8.0, K, = 4.0) can be obtained, while the heading deviation 1s 
increased a small amount. However, it appears that a rudder angle of + 32.6 degree 
and rudder speeds of + 16.3 deg/sec are required to achieve this. 

In practice there are limitations on both the maximum rudder angle and the 
maximum rudder speed. The maximum rudder angle hmits the maximum moment 
which can be applied. Figures 4.4 through 4.8 compare the behavior of the open 
system with that of the rudder roll stabilized ship when the rudder angle is limited to 
+ 20 degrees and the rudder rate limit is + 10 deg’sec. Each Figure shows the effect 
of a different compensation design. The svstem of Figure 4.5 appears to have 
minimum rudder activity and minimum yawing, but (from Table 5, Appendix I) the 


maximum roll angle is slightly larger than for the other design. 
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Data for various real rudder rate limits and various rudder angle limits were 
obtained and are given in Table 5. From this table, we can summarize the results as 


shown in Table 3: 


TABiIEiSs 


MINIMUM VALUES OF ® WITH DIFFERENT ZERO LOGA TGS 


NTIAX 


@ = 0.5 


een] 2 |e | [he [On 
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[300 | 10 | a0 [| 40 | 12 [90 
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It is clear that considerable reduction in roll angle can be achieved. A major 
requirement is high rudder rate. However for a given rudder rate good roll reduction 
can be obtained with several locations of the compensator zero. It may also be 
observed that roll reduction is accompanied by increased yawing unless the rudder rate 
is quite high. Thus a number of options are available. 

Maximum rudder rate is determined by the design of the rudder actuator system, 


but rudder rate less than the maximum could be set in the autopilot. In like manner 
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rudder angle limits can be set in the autopilot, and with proper design the feedback 
gains K, and K, can be made adjustable. Thus it 1s possible to design a rudder roll 
stabilization system which has considerable capability for roll reduction, but permits 
the user to choose the specific operating condition. 

Figures 4.9 through 4.12, shows system response for different feedback gains (K, 
and K,) with the same zero location when a complex wave causes rolling. Rudder 
activity and yaw motion are compared with that which would result if no stabilization 
Was used. 

Figures 4.13 through 4.16 show the effect of rudder rate on the maximum roll 
angle when the maximum rudder angle is limited. It 1s clear that for small rudder angle 
limits there is very little reduction in roll angle, and the rudder rate and zero locations 
have no noticeable effect. As the rudder angle limit is increased greater roll reduction 
becomes available providing the rudder rate is increased. From Figure 4.16 it is seen 
that maximum roll reduction is obtained when both rudder angle and rudder rate are 
maximum. Under these conditions the location of the compensator zero 1s important 
and from the data available it appears that the zero locations should be in the range of 
ie Ze 2.0. 

Figures 4.17 through 4.20 give a simulation study of the effects of roll 
stabilization on ship yawing. It is seen that those conditions which provide maximum 
roll reduction also result in increased yawing. Thus the user must choose operating 
values which best suit his needs in a particular sea state. 

Design of the roll stabilization was accomplished by considering rolling motions 
at the natural frequency only. To verify the results, the system was tested at a number 
of frequencies within the wave spectrum. Results are illustrated in Figures 4.21 
through 4.25. Rolling is maximum when the system is ‘open system’ ie, no 
Stabilization. Suppression of roll is greatest with an ideal rudder. When the rudder 
angle and rate are limited the effectiveness of the stabilization is reduced. A high 
rudder rate is essential, and if this is combined with a large rudder angle, roll reduction 


is Significant over the entire spectrum. 
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B. 


SIMULATION RESULTS 
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0.5 rad 


(t) 


SH 


foc. 


30d 


Cg 


RUDDER RATE (0BG/SEC) 


D:Z=2.0 


C:Z21.6 





i 3 
a2 2 : 
_ : 
2 
3 

qo 7 sd 80 fo 

(a0) JENV ONKIVGH 
a d 
d d 
J i 
Ee s ; 
die cs are 
H < _ 
g 3 
3 go #70 @0o #8 +3 





Figure 4.17 Maximum Rudder Rate vs Maximum Heading Angle, 


Maximum Rudder Angle = 5 deg 
0) = 0.5 rad/sec. 
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Figure 4.18 Maximum Rudder Rate vs Maximum Heading Angle, 
Maximum Rudder Angle = 10 deg 
a) = 0.5 rad/Scec. 
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Figure 4.19 Maximum Rudder Rate vs Maximu 
Maximum Rudder Angle = 20 d 
w = 0.5 rad'Sec. 
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Figure 4.20 Maximum Rudder Rate vs Maxim 
Maximum Rudder Angle = 30 
G) = 0.5 rad'sec. 
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Figure 4.21 Performance of Roll Reduction : Fr cepa vs Maximum Roll Angle 
Z=1.0, K 2.0, K 4.0 


Open Sy “stem | vs Ideal’Rudder. 
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Figure 4.22 Performance oy SO Reduction : Frequency vs Maximum Roll Angle 
Real Rudder : Maximutha Rudder Rate = + 5 deg/sec. 
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Figure 4.23 Performance of Roll Reduction : Frequency vs Maximum Roll Angle 


Real Rudder : Maxine Rudder Rate = + 10 deg/sec. 
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Figure 4.24 Performance ae Sea Say vs Maximum Roll Angle 
Real Rudder: Maximurh Rudder Rate = + 15 deg/sec. 
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Figure 4.25 Performance of a oe Frequency vs Maximum Roll Angle 


Real Rudder : ‘aad Rudder Rate = + 20 deg/sec. 
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V. CONCLUSION AND RECOMMENDATION 


CONCLUSION 


From the results of the Rudder Roll Stabilization (RRS) computer simulations, 


the following conclusions have been reached : 


(1) 


(2) 


(3} 


(4) 


(5) 


Simple models can be derived to describe how wave forcing input and rudder 
angle affect the roll angle. 

A properly designed roll controller will minimize the roll angle response and 
provide the desired heading angle as well. 

This thesis has presented some practical aspects of a Rudder Roll Stabilization 
system. It has been shown that the steering machine may act as a severe 
limitation. When the rudder speed 1s sufficiently increased good roll reduction 
can be realized on existing ships. When it is possible, by means of a modified 
rudder construction, to increase the effect of the rudder on the roll without 
changing its effect on the heading an even better performance can be 
obtained. 

It appears that both the rudder angle limit and the rudder speed limit have to 
be chosen carefully in order to realize a reasonable roll reduction. That is, the 
effectiveness of a RRS svstem is directly dependent upon available rudder 
moment. The available moment, in turn, depends on rudder rate. RRS 
system performance is rudder rate limited. Another parameter which 1s very 
important with respect to the roll reduction achieved by means of the rudder 1s 
the gain of roll and roll rate feedback (see Tables 3, 5 and 7). 

The primary objective of the RRS has been achieved. That is, maximum roll 
reduction up to 55 percent has been measured with an ideal rudder (Z= 1.0, 
Kk, =4.0, K,=4.0), and maximum roll has also been reduced to 47.5 percent 
with a reasonable real rudder (Z=1.0, K, =4.0, K,=4.0, 6, =10 deg/sec, 
Onax 7 20 deg), (see Tables 4 and 5). 
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B. 


RECOMMENDATION 


The results and conclusions of the RRS computer simulation lead to the 


following recommendation : 


(1) 


(2) 


(4) 


(5) 


The results of this thesis demonstrate that roll stabilization using the rudder is 
feasible. Only one roll controller has been studied. [t used roll angle and roll 
rate feedback. [It is recommended that roll controllers including roll 
acceleration feedback should be studied. 

In this thesis it has been assumed that cross coupling between the parameters 
of roll and yaw rate ( N@® ) is zero. However, it is recommended that cross 
coupling between these parameters be considered in further studies on this 
topic: 

The study does not consider wave input effects on the yawing moment. 
Further research should consider these effects. 

Future research should investigate the use of advanced adaptive control and 
optimal estimation techniques to solve the problem of Rudder Roll 
Stabilization. 

Future research should investigate the effects of various sea state conditions 
on the Rudder Roll Stabilization. 
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APPENDIN A 
CONCEPT BEHIND THE RUDDER-ROLL-STABILIZATION 


Roll stabilization is actually a secondary function of the rudder. Its primary 
purpose is to stecr. Hlowever, as shown in Figure A.I, the rudder produces 
simultaneously a roll moment and a yaw moment which 1s needed to change ship 
course. 

Typical ship response periods to a roll moment are about 8 to 12 seconds, where 
as, tvpical response periods a yaw moment are about 30 to 35 seconds. 

The significant diflerence between these response periods allows the simultaneous 
superposition of yaw and roll control signals on the rudder without adversely aflecting 
the response in either mode. 

The rudder moment is proportional to of the ship velocity and to the rudder 
angle. Using the rudder for roll stabilization involves opposing the wave induced roll 
moment by the rudder induced roll moment. 

Stabilization is achieved by adding to the steering control signal a roll control 


signal, 


Yaw ARM 


VAW MOMENT = LIFT x ARM 
YAW AXIS 


| 





Figure A.l Rudder Moments and Ship Response. 
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APRENDEX B 
THE ROLL MODEL WITHOUT COMPENSATOR 





Wf 


Figure B.1) Block Diagram Model of Ship’s Dynamics with Autopilot 
(without compensator). 
where 
G 
G 
G 


, 1s the integration function (= 1/S}. 


5 is the ship’s heading equation {= (St + 1)}. 


ed 


is a constant relating rudder angle to heading rate {= Ng}. 
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G, is the rudder servo equation {= KRSER/(S + KRSER)}. 
G, is the autopilot equation (= KHP + KHD*S}. 

G, is the ship’s roll equation (=o 7/(s? + 260 s+ )}. 

G, isa constant relating rudder angle to roll moment {= K¢}. 
Gg is a constant relating rudder angle to heading rate {= Ng}. 
Gy 1s a constant relating heading rate to roll moment {= K,}. 
C is a damping ratio of the roll dynamics. 

(is the natural frequency of the roll dynamics. 

KRSER — is a constant relating rudder servo gain. 

KHP is a constant relating heading gain (Autopilot). 
KHD is a constant relating heading rate gain (Autopilot). 


letting w,=0.0, 


For convenience, products of the types G,G,G,... are written as G 


easy) 


P=G,(m@, + A-B) 


A= - Gy 9457) 
B= Gy 9D 
OF 913 3,45? 
D=C=E 
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(B.1) 


(B.2) 


(B.3) 


(B.4) 


(B.3) 


Thus 
D= -G)4345D - Gg 
_ - Go® 
L+G) 4345 


D= GC) - Gi 5457) - Gy 9D ) 


(G Ge 
= G, @, a 1 .2,.4.5,.7,8 20 (:p 
1,2,3,4,5 


a Go(l+G, 33.45) 


MO, 1+G) 54345 - Gi 2456.78 - S2,6,8,9 


where Gg is zero 


(B.6) 


(B.7) 


(B.8) 


(B.9) 


(B.10) 


The physical interpretation of this is that in a beam sea, the waves do not cause 


yawing, and there is no cross coupling from roll to yaw, so the ship stays on course 


and the autopilot does not need to give rudder commands. Thus ship rolling is not 


affected by the rudder, because the rudder remains in the neutral position. 
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APPENDIX C 
THE ROLL MODEL WITH COMPENSATOR 





igure C.1 


Block Diagram. Model of Se Dynamics with Autopilot 
(With Compensator). 


where 
G 
G 
G 


, is the integration function {= 1/S}. 
, 1s the ship’s heading equation {= (St. + 1}. 


3 18 a constant relating rudder angle to heading rate {= Ng}. 
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G, is the rudder servo equation {= KRSER/(S+ KRSER)}. 

G, is the autopilot equation { = K BLP deeitetiee Si 

G, is the ship's roll equation {= 0 7/(s?-+ 26a, s+ *)}. 

G. isa constant relating rudder angle to roll moment (~Nae 
Gg isa constant relating rudder angle to heading rate L=Ng}- 
Gg 1s a constant relating heading rate to roll moment {= K}}. 
G, is a roll compensated equation {K,*S+ K,}. 

C is a damping ratio of the roll dynamics. 

@, is the natural frequency of the roll dynamics. 


KRSER © is a constant relating rudder servo gain. 
KHP is a constant relating heading gain (Autopilot). 
KHD is a constant relating heading rate gain (Autopilot). 


letting @, = 0.0, 


For convenience, products of the types G,G,G,... are written as Gy 9 30 


P=G,(@, + A-B) (ei 
A= -G)45,D-G,,G® (Ge 
Be GoD (C.3) 
C= -G,5345D-G,,G (C.4) 
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E= G.® (26) 
nus 
D= -G 


1,2,3,4,3? : G34G,@ ; Gg? 


Z aCe oat G.)® 





eu) 
I+G) 9345 

P= GC @, - Gy 4457) - Gy7G,® - G, 4D) 

Gr oO, ie (C1 2.3.4.4 570. G124578 + 22 3,490. * G28 9) 7 Gres G } @D (C.8) 

1+G Te 
1,2,3,4,5 
pee OS) eee (C.9) 
MO, 1G) 5345 - G1 2.45,6,7.8 - 2,689 7 GlG4.6.7-G23.4.6.9) 
where Ge is Zero 

el GAl+G, 2,3,4,5) (C.10) 
MO, I+Gi 5345 + GlGy67-G23469) 
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APPENDIX D 
COMPUTER PROGRAM FOR ROLL, UNCOMPENSATED 


TITLe Model Without Compensator 


KAKKKKKKKKKRKKKKKKKKKKKKRKKRKRKRKRRRKRRRKRRRRRRKRRRKRKRRKRKRRKRRRRKRRRRKRRRRRR 


AKX KKK 
KK The Roll Response of the Open System RAK 
kK KKK 


KAKKKEKKKKKKREKKKRKRREKRKKRRKRRRRRRRRKRRRRRKRKRKRRRRKRRRKRRRKRRERRRRRRRRRRKRKI 
* 


CONST ZETA=0.18, OMEGAN=0.59, 
KWH=0.0, KWR=1.0 


PARAM W=0.3 
PARAM A0=0.06 


* 
INITIAL 
YOQ= 


=0.0 
YDO=0.0 
* The Wave Equation 


TH=W*TIME 
WAVE=A0*SIN(TH) 


x The Roll Equation 


SO SAP ETA OnECLt ar? SUM 
SUM=(2*ZETA*OMEGAN*RORATE )+(OMEGAN**2*ROLL ) 
ROACC=ERROR 
RORATE= =INTGRL(YO, ROACC) 
roranq=rorate 180. 1S) ale! 
ROLL= NTGRL(YO, eae 
i ROLANG=ROLL 180. /3.14 


CONTRL FINTIM=500.0, DELT=0.1 
PRINT 1.0,rorang,ROLANG 
SAVE Q.1,rorang,ROLANG 


END 
PARAM W=0.4 
PARAM A0=0.10 


END 
PARAM W=0.5 
PARAM A0Q=0.15 


END 
PARAM W=0.6 
PARAM A0=0.10 


END 
PARAM W=0.7 
PARAM A0=0.09 


END 
PARAM W 
PARAM A 


END 


PARAM W=0.9 
PARAM AO=0.07 
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* 


END 
PARAM W=1.0 
PARAM A0=0.05 


END 
PARAM W=1.1 
PARAM A0=0.04 


* 

GRAPH iG DE=TEK618,P0=0, NI=8) TIME(UN=SEC) ROLANG(UN=DEG) 
LABEL (G) ROLL ANGLE 

END 

STOP 


71 


APPENDIX E 


COMPUTER PROGRAM FOR ROLL COMPENSATED, IDEAL RUDDER 


TITLE Model Ideal Rudder With Compensator : Sine Wave and Complex Wave 


KEKKKKKKKKKKKKRKRKEKKRKKRRKKRERKKRERKKKKRRRKRKKKKRKRKRRKRKRKKKRRRKRRKRKKKRRKK 


kkk kkk 
KKK The Roll Response of the Open System and ARK 
k*k* kK 
KKK Ideal Rudder with Roll and ARK 
kkk kk 
KK Roll Rate Feedback AKK 
kk * kk 


KKEKAKKKKKKKKKKKKKKEKKKKKKAKKKRKAKKRKKKRKKKRKKKKRKRKKRKKRKKKKRKRKKKRKKRKRKKKKK 
* 


TITLE THE ROLL RESPONSE WITH COMPENSATOR : IDEAL RUDDER 


CONST COMMAX=1000.0, COMMIN=-1000.0, RATMAX=1000.0, ee 1000.0 aa 


KPSDR=0.077, KPHDR=0.20, KPHSDT=4.9, KPHPH=0.0, 
KWH=0.0, KWR=1.0,... 

KRSER=10.0, TPS= 6. 0y 

OMEGAN=0.59, ZETA=0. 18, 

Z=0.025, P= 1000, 

KHP=40000.0, RHD=0. iO). 

KRP=6.0, KRD=6.0, KRDD=0 .0 


W1=0.3, W2=0.4, W3=0.5, W4=0.6,... 
W5=0.7, W6=0.8, W7=0.9, W8=1.0,... 
A1l=0.06, A2=0.1, A3=0.15, A4=0.1, 
A5=0.09, A6=0.07, A7=0. 07, A8=0. 05, 
K=0.3 

k 

INITIAL 
YO=0.0 
YDO=0.0 

DERIVATIVE 
HEDCOM=0 .0 


ERROR1=HEDCOM-YAW 
x The Auto-Pilot Equation 


COMK=ERROR1*KHP 
COM=ZEROPL(Y0,Z,P,COMK) 
ERRORZ=COM=GC 
COMLIM=LIMIT (COMMIN , COMMAX, ERRORZ) 
ERROR3=COMLIM-DELTA 


* The Rudder Servo Equation 


ee ee 
DELTA=INTGRL( YO, DELLIM) 

RUDRAT=DELLIM 

RDRANG=RUDRAT*180./3.14 

RUDDER=DELTA 

RUDANG=RUDDER*180./3.14 


The Wave Equation 
PURE SINE WAVE 


TH=W3*TIME 
WAVEK=A3*SIN(TH) 


HO OO OO 


ue 


‘i COMPLEX WAVE 


WAVE=A1*SIN ere tnexcen exrness. 
A4*SIN(W44TIME )+A5*SIN(WSATIME )+A6*SIN(WO*TIME) +... 
AT*SIN(W7*TIME ) +AS*SIN (WE*TINE ) 

WAVEK=WAVE*K 


* The Ship's Yaw Equation 


SUM1=DELTA*KPSDR-ROLL2*KPHPH+WAVEK*KWH 
YAWR=REALPL (YO ,TPS,SUM1) 
YAW=INTGRL (YO, YAWRS 
YWRANG=YAWR*180./3.14 
YAWANG=YAW*180./3.14 


THE SHIP'S ROLL EQUATION (OPEN SYSTEM) 


SUMO=( 2*ZETA*OMEGAN*RORATI ) + (OMEGAN**2*ROLL1 ) 
ERRORO= (OMEGAN**2*WAVEK*KWR ) - SUMO 
ROACC1=ERRORO 

RORAT1=INTGRL(YO,ROACC1) 

ROLLI=INTGRL (YO, RORAT1 ) 
ROLAG1=ROLL1%*180./3.14 


THE SHIP'S ROLL EQUATION (COMPENSATED SYSTEM) 


SUM2=DELTA*KPHDR- YAWR*KPHSDT+WAVEK*KWR 
See eee nO DNAR Ree )t (OMEGAN**2*ROLL2 ) 
ERROR4= (OMEGAN** 2*SUMZ ) -SUM3 

ROACC2=ERROR4 

RORAT2=INTGRL( YO, ROACC2) 
ROLL2=INTGRL (YO ,RORATZ) 
RORAG2=RORAT2*180./3.14 
ROLAG2=ROLL2*180./3.14 


The Roll Feedback Equation 
e GC=KRP*ROLL2+KRD*RORAT2 


TERMINAL 

NELnOD STIFF 

CONTROL FINTIM=400, DELT=0.1 

PRINT 1.0, RDRANG, RUDANG, YAWANG ,ROLAGI , ROLAGZ 
SAVE Q.1, RDRANG, RUDANG, YAWANG,ROLAG1 ,ROLAGZ 


GRAPH (TOP/G,DE=TEK618,PO=0.0,6.5, NI=6) TIME(UN=SEC) RUDANG(LE=1.5,... 
UN=DEG , LO=-30 , SC=10 ,NI=6.0 

LABEL a ‘IDEAL SYSTEM' SINE WAVE INPUT 
LABEL (TOP K1=4.0, K2=4.0, W=0.5 RAD/SEC 


GRAPH Meee nee ecb TERGTS (FO=0 1000, NI=6) TIME(UN=SEC) YAWANG... 
(LE=1.5,UN=DEG ,LO=-3.0,SC=1.0,NI=6.0) 

LABEL pre ‘IDEAL SYSTEM' SINE WAVE INPUT 

LABEL MIDDLE K1=4.0, K2=4.0, W=0.5 RAD/SEC 


GRAPH (BOTTOM/G,OV,DE=TEK618,P0=0.0,1.0, NI=5) TIME(UN=SEC) ROLAG1... 


+ + + 


+ + 


+ FOF 


(LO=-25 ,SC=10.0,NI=5,LE=1.5,UN=DEG) ROLAG2(LO=-25 ,SC=10.0,NI=5,LE=1.5,... 


UN=DEG 

LABEL tees 'OPEN SYSTEM AND IDEAL SYSTEM' SINE WAVE INPUT 
LABEL BOTTOM K1=4.0, K2=4.0, W=0.5 RAD/SEC 

END 

STOP 


ie 


APPENDIX F 
COMPUTER PROGRAM FOR ROLL COMPENSATED, REAL RUDDER 


TITLE Model Real Rudder With Compensator : Sine Wave and Complex Wave 


KAKKKKKKKKKKKKKAKAKKKKKKRKKKEKRKRKKKKRRKKKKKKKRRKKRKEKKKKRKKKKKKEKKKKKKKRKK 


KKK Lk 
AK The Roll Response of the Open System and KKK 
KK te 
na Real Rudder with Roll and Roll Rate KKK 
ka KK 
ARK Feedback, and Nonlinearity in the AK 
KK dhe 
gers Rudder Servo. an 
KKK Kk 


KKEKKKKKKKKKKKKRKKKEKKRKKKKKRKKRKRKKKKKKRKERKKRKKKKKRKRKEKKKKKRKKKRKKRKRKKKK 


TITLE THE ROLL RESPONSE WITHOUT/WITH COMPENSATOR "K1=4.0, K2=4.0"! 
x OPEN SYSTEM AND GOMPENSATED = Sis fen 
‘A IDEAL RUDDER 
XONST COMMAX=0.0873, COMMIN=-0.0873, RATMAX=0.1745, RATMIN=-0.1745,... 
XONST COMMAX=0.1745, COMMIN=-0.1745, RATMAX=0.1745, RATMIN=-0.1745,.. 
XONST COMMAX=0.3491, COMMIN=-0.3491, RATMAKX=0.1745, RATMIN=-0.1745,.. 
XONST COMMAX=0.5236, COMMIN=-0.5236, RATMAX=0.1745, RATMIN=-0.1745,.. 
XONST COMMAX=0.0873, COMMIN=-0.0873, RATMAX=0.2618, RATMIN=-0.2618,.. 
XONST COMMAX=0.1745, COMMIN=-0.1745, RATMAX=0.2618, RATMIN=-0.2618,... 
XONST COMMAX=0.3491, COMMIN=-0.3491, RATMAX=0.2618, RATMIN=-0.2618,... 
XONST COMMAX=0.5236, COMMIN=-0.5236, RATMAX=0.2618, RATMIN=-0.2618,... 
XONST COMMAX=0.0873, COMMIN=-0.0873, RATMAX=0.3491, RATMIN=-0.3491, 
XONST COMMAX=0.1745, COMMIN=-0.1745, RATMAX=0.3491, RATMIN=-0.3491, 
XONST COMMAX=0.3491, COMMIN=-0.3491, RATMAX=0.3491, RATMIN=-0.3491, 
XONST COMMAX=0.5236, COMMIN=-0.5236, RATMAX=0.3491, RATMIN=-0.3491, 
XONST COMMAX=0.0873, COMMIN=-0.0373, RATMAX=0.0873, RATMIN=-0.0873,... 
XONST COMMAX=0.1745, COMMIN=-0.1745, RATMAX=0.0873, RATMIN=-0.0873,... 
XONST COMMAX=0.3491, COMMIN=-0.3491, RATMAX=0.0873, RATMIN=-0.0873,... 
CONST COMMAXK=0.5236, COM eo) eaeer RATMAXK=0 .0873, a -0.08 7370 

KPSDR=0.077, KPHDR= 20, KPHSDT= 4.9, KPHPH=0.0, 

KWH=0.0, KWR=1. Ousenne 

KRSER= 10. 07 SEPS= 6. Ti 

OMEGAN=0. 59, ZETA=0. 18, 

Z=0.025, P= 1000, 

KHP=40000. OF KHD= 0. 0, 

W1=0.3, W2= 0. 4, W3=0. 5, “W4=0.6, 

W5=0.7, W6=0.8, W/=0.9, W8=1.0,... 

Al=0.06, A2=0.1, A3=0. 15, A4=0.1,. 

oe Ao=0. 07, A7=0. 07, A8=0. 05,. 
- KRP=4.0, “KRD= 4.0, KRDD=0.0 
INITIAL 

YO=0.0 

YDO=0.0 
DERIVATIVE 

HEDCOM=0 .0 
" ERROR1=HEDCOM-YAW 
A The Auto-Pilot Equation 


COMK=ERROR1*KHP 
COM=ZEROPL(Y0O,Z,P,COMK) 
ERROR2=COM-GC 
COMLIM=LIMIT(COMMIN , COMMAX , ERRORZ) 


714 


ERROR3=COMLIM~-DELTA 
* The Rudder Servo Equation 
RPE TN BA DMAe KRSER“ERRORS ) 
DELTA=INTGRL(YO ,DELLIM) 
RUDRAT=DELLIM 
RDRANG=RUDRAT*180./3.14 
RUDDER=DELTA 
RUDANG=RUDDER*180./3.14 
The Wave Equation 
PURE SINE WAVE 


TH=W3*TIME 
WAVEK=A3*SIN(TH) 


COMPLEX WAVE 
pene eect snexoin neteiwe st. 


A Oe A A Oe OO OO Ot 


A4X*SIN(W44TIME )+A5*SIN(WS*TIME )+AG*SIN(WO*TIME)+... 
AT*SIN(W7*TIME ) +A8*S IN (WS*TIME ) 
WAVEK=WAVE*K 


* The Ship's Yaw Equation 


SUM1=DELTA*KPSDR~ROLL2*KPHPH+WAVEK*KWH 
YAWR=REALPL(YO,TPS,SUM1) 
YAW=INTGRL(YO, YAWR 
YWRANG=YAWR*180./3.14 
YAWANG=YAW*180./3.14 


* 

* THE SHIP'S ROLL EQUATION (OPEN SYSTEM) 
SUMO=( 2*ZETA*OMEGAN*RORAT1 )+(OMEGAN**2*ROLL1) 
ERRORO=(OMEGAN**2*WAVEK*KWR) -SUMO 
ROACC1=ERRORO 
RORAT1=INTGRL (YO, ROACC1) 
ROLL1=INTGRL (YO, RORAT1 ) 

, ROLAG1=ROLL1*180./3.14 

* THE SHIP'S ROLL EQUATION (COMPENSATED SYSTEM) 


SUM2=DELTA*KPHDR~ YAWR*KPHSDT+WAVEK*KWR 
Se ey gee Oa RR CRAT? |. COMBGAN’*2*ROLL2 ) 
ERROR4= (OMEGAN**2*SUM2 ) ~SUM3 
ROACC2=ERROR4 
RORAT2=INTGRL( YO ,ROACC2) 
ROLL2=INTGRL(Y0O,RORAT2) 
RORAG2=RORAT2*180./3.14 

a ROLAG2=ROLL2*180./3.14 


‘i The Roll Feedback Equation 
F GC=KRP*ROLL2+KRD*RORAT 2 


TERMINAL 

WITIG (CVE) SILI gy 

CONTRL FINTIM=400.0, DELT=0.1 

PRINT 1.0,RDRANG, RUDANG, YAWANG, ROLAG1 , ROLAG2 
SAVE 0.1,RDRANG,RUDANG, YAWANG, ROLAGI , ROLAG2 


GRAPH (TOP/G,DE=TEK618,P0=0.0,6.5, NI=6) TIME(UN=SEC) RUDANG(LE=1.5,... 


UN=DEG , LO=-30,5C=10,NI=6.0 
LABEL op 'REAL SYSTEM! SINE WAVE INPUT 
LABEL ine P K1=4.0, K2=4.0 


GRAPH (MIDDLE/G,OV,DE=TEK618,P0=0.0,4.0, NI=6) TIME(UN=SEC) YAWANG... 


Js 


(LE=1.5,UN=DEG, LO=-3.0,SC=1.0,NI=6.0) 
LABEL we ‘REAL SYSTEM! SINE WAVE INPUT 
LABEL (MIDDLE K1=4.0, K2=4.0 

| 


GRAPH (BOTTOM/G,OV,DE=TEK618,P0=0.0,1.0, NI=5) TIME(UN=SEC) ROLAG1... 
(LO=-25 ,SC=10.0,NI=5,LE=1.5,UN=DEG) ROLAG2(LO=-25,SC=10.0,NI=5,LE=1.5,. 
UN=DEG 

LABEL SOTTO ‘OPEN SYSTEM AND REAL SYSTEM' SINE WAVE INPUT 

LABEL BOTTOM K1=4.0, K2=4.0 

END 

STOR 
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APPENDIX G 
COMPUTER PROGRAM FOR PLOTTING IN DISSPLA: A 


KARKAKKAKRAKKRKKKKKKRKRKKKKKKKKKRKKKRARKKKKKKRRKAKRKKKKK 


Ka KK 
ARK This is a program for plotting ARK 
KX KK 
KK in DISSPLA in maximum rudder rate ARK 
KAK sock 
KKK VS maximum roll angle and maximum ARK 
KK KK 
ARK rudder rate vs maximum heading angle *** 
KA toe 


KRREEKKKKKKKKKKKRKRKKRKKRKRKKKRKKRKAKRKAKKKKKKRKKKKKKAKAKKKKKK 


Cc 
C PLOTTING THE GIVEN FUNCTION 
PROGRAM GRF2D 
CALL COMPRS 
CALL XYPLOT 
CALL DONEPL 
SLOP 
END 


SUBROUTINE XYPLOT 
DIMENSION X(50),Y(50) 
REAL X,Y 

INTEGER I,J,N,M 

CALL PAGE(21.,18.) 


c 
EeGRE ist 


CALL PHYSOR 2. 12, :) 

CALL AREA2D(7.,5 

CALL YAXANG 20. it 

CALL HEIGHT .25) 

CALL COMP 

CALL ENANE ‘RUDDER RATE (DEG/SEC) | ,21) 
CAML ae ‘HEADING ANGLE (DEG)' 19) 


5 

4 

20 CONTINUE 
CALL CURVE(X,Y 
CALL RESET iy 
CALL ENDGR(0) 


C 
C GRF 2ND 
Cc 


CAEL PawSOR(1Z: ,12.) 

CALL AREAZD(7.,5. 

CALL ae) RUDDER RATE (DEG/SEC) ' aay 
CALL YNAME('HEADING ANGLE (DEG)' ,19) 
CALL DASH 

CALL GRAF me i" Noe Om, 0. lO). 5 ) 

CALL GRID 

CALL RESET | SASH") 

aie icy ie . ‘Oy? 


i 


READ(5,25)X(J) ,¥(J) 
25 FORMAT(F4.1, F6.1) 
30 CONTINUE 
CALL CURVE(X,Y, 4 0) 


CALL Re oe 
CALL ENDGR(O) 


CALL See one 
CALL AREAZD(7.,5. 

CALL set ete RATE EE ae ey 
CALL YNAME('HEADING ANGLE (DEG)',19) 
CALL DASH 


READ(5, 35 

35 FORMAT (F4. 
40 CONTINUE 
CALL CURVE(X,Y 

CALL RESET('TH 

CALL ENDGR(0) 


C 
C GRF 4TH 
Cc 
CALL BAER 2D (7< 5.5) 





CALL AREA2D(7.,5.) 
G CALL HEIGHT(0.2) 
CALL XNAME('RUDDER RATE (DEG/SEC)! ce 
CALL YNAME('HEADING ANGLE (DEG)',19) 
C CALL RESET('HEIGHT' ) 
CALL DASH 
GAIT rs Do) 2D Oe, 0 ole Ome 





CALL GRrDt lel 
C CALL HEIGHT (0. 05) 

CALL RESET('DASH' ) 

CALL THKCRV(0O. 05) 

DO 50 J =1, 
READ(S, a8)x(0),¥(9) 

45 FORMAT(F4.1, F6.1) 
50 CONTINUE 

CALL CURVE(X,Y,4,0) 

CALL RESHE UAEL! } 

CALL ENDGR(0O) 


MESSAGE 


CALL PHYSOR 
CALL AREAZ2D 
CALL HE@Gue 
CALL COMPLX , 





aA 


SNS 





5 


,13,4 0) Or 3238) 


] 
Z ales 14. my 2 
] 13 


? / 


B 

C 

D 13s iy , 

A; MAXIMUM RUDDER ANGLE : 5 (DEG) 
B; MAXIMUM RUDDER ANGLE : 10 bee 
D 
ie 


Wo 
NRF OO 
OMmO vi 


IN 


MAXIMUM RUDDER ANGLE : 20 (DEG 
NAXIMUM RUDDER ANGLE : 30 (DEG 


AA ME 
‘2 
Ned 
c 
cS 
a 
OC] 
wm 
Ww 
ed 
© 
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APPENDIX H 
COMPUTER PROGRAM FOR PLOTTING IN DISSPLA : B 


KRAEKKKKKKKKKKKKRKRKKRKRKRKRKKRKRKRRKKRKRKRKKKKKAKRKKKRKRKAKKRRRKRRER 


AK KA 
KA te 3 KK 
aA This 1s a program for plotting ee 
ARK in DISSPLA in FREQUENCY VS RAK 
KK KK 
ARK maximum roll angle OF OPEN SYSTEM, RAK 
KKK KAK 
ARK IDEAL RUDDER AND REAL RUDDER. ARK 
KAK sok ® 


KRAEKRKEKKKKRKRKAKRKKRKRKRKRKRKRKKKARKKRRRRRKRKARKAKKRARKRRRRRKKKKKRKREKA 


C 
ec GRAPH OF AY DATA GIVEN 
PROGRAM GRF2ZD 
CALL TEK618 
CA en © LOT 
CALL DONEPL 
Sar 


END 
C PLOTTING THE GIVEN FUNCTION 
SUBROUTINE XYPLOT 
DIMENSION N x(100), *y1(100), Y2(100) , IPKRAY(100) 
REAL X,Y 
eee 7 "IN M 
C SET PAGE AND SUBPLOT SIZES 
CALL PAGE(15.,11. 
CALL AREA D(? oo 
CALL HEIGHT 25) 
C LABEL AND DRAW AXIS 
CALL NAME ( (RAD/SEC)' 11) 
CALL YNAME('ROLL-MAX (DEG)' 14) 
CALL COMPLX 
CALL FRAME 
CALL GRAF (0.2 
CALL GRID(1,15 
DO 20 I 


CALL SPLINE 
READ(S, S8)x(x), YL(I), Y2(I) 
55 FORMAT (F3.1,F7.1,F6.1) 
20 CONTINUE 
C SET UP FOR LEGEND 
MAXLIN = LINEST(IPKRAY,100,15) 
CALL LINES ee SYSTEMS! , IPKRAY, 1) 
CALL LINES ('IDEAL SYSTEMS! , IPKRAY, 2) 
CALL LEGLIN 
C PLOT CURVES 
CALL CURVE(X, Y1,9,1) 


Os. 200.0440, 22.10) 


CALL DASH 
CALL CURVE (X /¥2,9,1) 
CALL RESET('ALL') 


C WRITE LEGEND AND TITLE 
Gane coo EO CERES 5) 
CALL HEADIN('OPEN SYSTEM VS IDEAL SYSTEM $',100,-4,2) 
CALL ENDPL(0) 
RETURN 
END 


dg 


APPENDIX I 
SIMULATION RESULTS 


TABLE 4 
IDEAL RUDDER WITH DIFFEREMIRZ: tO Eee. ie> 





Tio [ao | ao | as [nr pa | 90 
30 [so | 40 [163 [6 | na | 8 














VAIS te es 
REAL RUDDER WITH DIFFERENT ZERO LOCATIONS 


Mm = 0.5 


Maximum Rudder Angle Limit = +5 degree 


Maximum Rudder Rate Limit = +5 degree/sec 


os [20 [40 [30 | 30 | 02 | 179 
[20 [so [40 | so [so] 03 | 180 













SO 


(aD > 
REAL RUDDER WITH DIFFERENT ZERO LOCATIONS (CONT’D.) 


Maximum Rudder Angle Limit = +10 degree 


Maximum Rudder Rate Limit = +5 degree/sec 


10 | a0 | 40 | 30 | 100 | 05 | 162 
[20 | so | 40 | so [100 [os 172 


Maximum Rudder Angle Limit = + 20 degree 










Maximum Rudder Rate Limit = +5 degree/sec 


Tos [20 | a0 | so | 32 [06 [iss 
Tis [60 | 40 | so | x6 [07 | 186 
[20 [so | 40 | so | ss [os [95 


Maximum Rudder Angle Limit = +30 degree 











Maximum Rudder Rate Limit = +5 degree/sec 


[os [20 | 40 | so | is2 [08 | 15s 
[20 [so | ao [so [mi | 09 | 93 














8 | 


ABILE.o 
REAL RUDDER WITH DIFFERENT ZERO LOCATIONS (CONT D.) 


Maximum Rudder Angle Limit = +5 degree 
Maximum Rudder Rate Limit = + 10 degree/sec 


i 
Tas [as [2s | ioo | 30) 03 [an 
a 
1s [60 [a0 [100 [so] 0s 
[20 [so | 40 | wo [so [03 | a7 


Maximum Rudder Angle Limit = + 10 degree 











Maximum Rudder Rate Limit = + 10 degree/sec 


i 
cis [as [23 [100 | i00 [06 | 143 
[20 [so | a0 [0 | ino] 06 [asa 


Maximum Rudder Angle Limit = +20 degree 












Maximum Rudder Rate Limit = + 10 degree/sec 


os [20 | 40 | 00 | 200 | 09 | aaa 
io [40 | a0 | 100 | 200 | 10 | 105 
20 [so | ao [100 [200 [0 | 7 
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TABLE 5 
REAL RUDDER WITH DIFFERENT ZERO LOCATIONS (CONT'D.) 


Maximum Rudder Angle Limit = +30 degree 
Maximum Rudder Rate Limit = +10 degree/sec 


os [20 | a0 | 00 | 3s9 [| a2 
Tis [60 | 40 | 100 [300 | 12 | 10s 
[20 [so [40 | 0 | 300] 12 | a8 


Maximum Rudder Angle Limit = +5 degree 















Maximum Rudder Rate Limit = +15 degree/sec 


“os [20 | 40 [150 [30 [03 | 181 


Maximum Rudder Angle Limit = +10 degree 









Maximum Rudder Rate Limit = +15 degree/sec 


ee ge as ee] ee | ase 
T10 [40 | a0 [iso | too [06 | as 
Tas [60 | 40 | 150 | 0 | 06 | 148 
[20 [so | 40 [150 | wo [06 [14s 
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TAMBIEE 5 
REAL RUDDER WITH DIFFERENT ZERO LOCA TIONS (COMED, 


Maximum Rudder Angle Limit = +20 degree 


Maximum Rudder Rate Limit = +15 degree/sec 


a 
20 [so [40 | 50 [0 [141 | 94 


Maximum Rudder Angle Limit = +30 degree 













Viaximum Rudder Rate Limit = +15 degree/sec 


A 
Tio [a0 [aos Daa | 50 


Maximum Rudder Angle Limit = +5 degree 












Maximum Rudder Rate Limit = +20 degree/sec 


[20 [so | 40 [| 20 | so [03 [at 
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MEAIEE SD 


REAL RUDDER WITH DIFFERENT ZERO LOCATIONS (CONT’D.) 


Maximum Rudder Angle Limit = +10 degree 


Maximum Rudder Rate Limit = +20 degree/sec 


Se eee mee 
os [20 | 40 | 200 | 00 | 06 | 160 
10 [40 | a0 | 200 | 100 [06 | 150 
Tas [60 | a0 | 200 | 100 [06 [a7 
20 [80 | ao [20 | 0 | 06 [iar 


Maximum Rudder Angle Limit = + 20 degree 
















Maximum Rudder Rate Limit = +20 degree/sec 


ST 
[60 | 40 | 20 | 200 [41 [93 
Oe 


Maximum Rudder Angle Limit = +30 degree 













Maximum Rudder Rate Limit = + 20 degree/sec 


Tos [20 | a0 [a | ma [i | 6 

A 
asso a0] 153] 00 [ne [0 
[20 [so | 40 [171 | 300 [20 | 59 
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TABLE 6 
IDEAL RUDDER WITH DIFFERENT GAIN 


Z 1.0 
@ 0.5 








20 [20 [99 | ws] 09 | 29 
Tso [so [am [36 | 16 [5s 





TABLE 7 
REAL RUDDER WITH DIFFERENT GAIN 


Z 
(D 


1.0 
0.5 
Maximum Rudder Angle Limit = +5 degree 


Maximum Rudder Rate Limit = +5 degree/sec 


so [so | so | so [os | 179” 
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Powe | 
REAL RUDDER WITH DIFFERENT GAIN (CONT'D.) 


Maximum Rudder Angle Limit = +10 degree 


Maximum Rudder Rate Limit = +5 degree/sec 


[so | so | so] oo [0s | 165” 


Maximum Rudder Angle Limit = +20 degree 











Maximum Rudder Rate Limit = +5 degree/sec 


[20 [20 | so | 3s | 06 | 161 
rso [so | so | isa [09 | 183” 


Maximum Rudder Angle Limit = + 30 degree 







Maximum Rudder Rate Limit = +5 degree/sec 


[20 [20 | so [iss [06 | 161” 


60 | 60 | so | 160 | 0s | 181 
[so | so | so | ise [09 [187 
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ABLE 7 
REAL RUDDER WITH DIFFERENT GAIN (CONT D.) 


Maximum Rudder Angle Limit = +5 degree 


Maximum Rudder Rate Limit = +10 degree/sec 


Maximum Rudder Angle Limit = +10 degree 














Maximum Rudder Rate Limit = +10 degree/sec 





Maximum Rudder Angle Limit = = 20 degree 


Maximum Rudder Rate Limit = +10 degree/sec 


“20 [20 | 99 | 9s [os | 29 
To [60 | 0 | 20 | 09 | 108 
Tso [so | too | 200 [10 | 108 











58 


WABLE 7 


hea DDER WITheDIFFERENT GAIN (CONT'D.) 


Maximum Rudder Angle Limit = + 30 degree 


Maximum Rudder Rate Limit = +10 degree/sec 


[20 [20 | 99 | 198 | 09 | 29” 
so | so | wo [0 | 12 | 108” 


Maximum Rudder Angle Limit = +5 degree 







Maximum Rudder Rate Limit = +15 degree/sec 


eo [oe 
A 
so | 60 | 150 | so] 03 [177 
so [80 | 150 [so [03 | an7 


Maximum Rudder Angle Limit = £10 degree 












Maximum Rudder Rate Limit = £15 degree/sec 


Tso [60 | 150 | 100 | 08 | 14s” 
[so [so | 150 [100 | os | 14s” 
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eel. 7 


REAL RUDDER WITH DIFFEREN Mente D) 


Maximum Rudder Angle Limit = + 20 degree 


Maximum Rudder Rate Limit = +15 degree/sec 


[so [so [iso | 20 | a | 935 


Maximum Rudder Angle Limit = + 30 degree 















Maximum Rudder Rate Limit = +15 degree/sec 






Maximum Rudder Angle Limit = +5 degree 


Maximum Rudder Rate Limit = + 20 degree/sec 


[so | 60 | 20 | so | 03 [177 
Tso | so | 200 | so | 03 | 177 
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TABLE 7 


REAL RUDDER WITH DIFFERENT GAIN (CONT’D.) 


Maximum Rudder Angle Limit = +10 degree 
Maximum Rudder Rate Limit = +20 degree/sec 


ae er 
[60 | 60 | 200 | 100 | 06 | 17 
[so [so | 200 | 100 | 06 | 147 


Maximum Rudder Angle Limit = + 20 degree 













Maximum Rudder Rate Limit = +20 degree/sec 


i a 
[20 [20 [99 | ws [0s | 29 
60 | 69 | 20 | 0 [1a | 95 
so [so | 20 [mo fair | 32 


Maximum Rudder Angle Limit = + 30 degree 








Maximum Rudder Rate Limit = +20 degree/sec 
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Rudder roll stabilization. 





